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We address the problem of spin-resolved scattering through spin-orbit nanostructures in graphene, 
i.e., regions of inhomogeneous spin-orbit coupling on the nanometer scale. We discuss the phe- 
nomenon of spin-double refraction and its consequences on the spin polarization. Specifically, we 
study the transmission properties of a single and a double interface between a normal region and 
a region with finite spin-orbit coupling, and analyze the polarization properties of these systems. 
Moreover, for the case of a single interface, we determine the spectrum of edge states localized at 
the boundary between the two regions and study their properties. 
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I. INTRODUCTION 

Graphene 1,2 — a single layer of carbon atoms arranged 
in a honeycomb lattice — has attracted huge atten- 
tion in the physics community because of many unusual 
electronic, thermal and nanomechanical properties. 3,4 In 
graphene the Fermi surface, at the charge neutrality 
point, reduces to two isolated points, the two inequiv- 
alent corners K and K' of the hexagonal Brillouin zone 
of the honeycomb lattice. In their vicinity the charge 
carriers form a gas of chiral massless quasiparticles with 
a characteristic conical spectrum. The low-energy dy- 
namics is governed by the Dirac-Weyl (DW) equation 5,6 
in which the role of speed of light is played by the elec- 
tron Fermi velocity. The chiral nature of the quasipar- 
ticles and their linear spectrum lead to remarkable con- 
sequences for a variety of electronic properties as weak 
localization, shot noise, Andreev reflection, and many 
others. Also the behavior in a perpendicular magnetic 
field discloses new physics. Graphene exhibits a zero- 
energy Landau level, whose existence gives rise to an un- 
conventional half-integer quantum Hall effect, one of the 
peculiar hallmarks of the DW physics. 

Driven by the prospects of using this material in spin- 
tronic applications, 7,8 the study of spin transport is one 
of the most active field in graphene research. 9-14 Sev- 
eral experiments have recently demonstrated spin injec- 
tion, spin- valve effect, and spin-coherent transport in 
graphene, with spin relaxation length of the order of 
few micrometers. 10,14 In this context a crucial role is 
played by the spin-orbit interaction. In graphene symme- 
tries allow for two kinds of spin-orbit coupling (SOC). 15 
The intrinsic SOC originates from carbon intra-atomic 
SOC. It opens a gap in the energy spectrum and con- 
verts graphene into a topological insulator with a quan- 
tized spin-Hall effect. 15 This term has been estimated 
to be rather weak in clean flat graphene. 16-19 The ex- 
trinsic Rashba-like SOC originates instead from inter- 
actions with the substrate, presence of a perpendic- 
ular external electric field, or curvature of graphene 
membrane. 16-18,20 This term is believed to be responsible 



for spin polarization 21 and spin relaxation 22,23 physics in 
graphene. Optical-conductivity measurements could pro- 
vide a way to determine the respective strength of both 
SOCs. 24 

In this article we address the problem of ballistic spin- 
dependent scattering in the presence of inhomogeneous 
spin-orbit couplings. Our main motivation stems from 
a recent experiment that reported a large enhancement 
of Rashba SOC splitting in single-layer graphene grown 
on Ni(lll) intercalated with a Au monolayer. 25 Further 
experimental results show that the intercalation of Au 
atoms between graphene and the Ni substrate is essen- 
tial in order to observe sizable Rashba effect. 26,27 The 
preparation technique of Ref. 25 seems to provide a sys- 
tem with properties very close to those of freestanding 
graphene in spite of the fact that graphene is grown on 
a solid substrate. The presence of the substrate does not 
seem to fundamentally alter the electronic properties ob- 
served in suspended systems, i.e., the existence of Dirac 
points at the Fermi energy and the gapless conical dis- 
persion in their vicinity. 

These results suggest that a certain degree of control 
on the SOC can be achieved by appropriate substrate 
engineering, with variations of the SOC strength on sub- 
micrometer scales, without spoiling the relativistic gap- 
less nature of quasiparticles. This could pave the way for 
the realization of spin-optics devices for spin filtration 
and spin control for DW fermions in graphene. An opti- 
mal design would require a detailed understanding of the 
spin-resolved ballistic scattering through such spin-orbit 
nanostructures, which is the aim of this paper. 

The problem of spin transport through nanostructures 
with inhomogeneous SOC has already been thoroughly 
studied in the case of two-dimensional electron gas in 
semiconductor heterostructures with Rashba SOC. 28-30 
Here the Rashba SOC 31 — arising from the inversion 
asymmetry of the confinement potential — couples the 
electron momentum to the spin degree of freedom and 
thereby lifts the spin degeneracy. In this case, a region 
with finite SOC between two normal regions has prop- 
erties similar to biaxial crystals: an electron wave inci- 
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Figure 1: (Color online) Illustration of the kinematics of the 
scattering at a N-SO interface in graphene. The circles rep- 
resent constant energy contours. 



dent from the normal region splits at the interface and 
the two resulting waves propagate in the SO region with 
different Fermi velocities and momenta. 28 This effect — 
analogous to the optical double-refraction — produces an 
interference pattern when the electron waves emerge in 
the second normal region. Moreover, electrons that are 
injected in an spin unpolarized state emerge from the SO 
region in a partially polarized state. 

Here we shall focus on the two simplest examples of SO 
nanostructures in graphene: (i) a single interface between 
two regions with different strengths of SOC; (ii) a SOC 
barrier, or double interface, i.e., a region of finite SOC 
in between two regions with vanishing SOC. 

Our analysis shows — in analogy to the case of 2DEG 
- that the ballistic propagation of carriers is governed 
by the spin-double refraction. We find that the scattering 
properties of the structure strongly depend on the injec- 
tion angle. As a consequence, an initially unpolarized 
DW quasiparticlc emerges from the SOC barrier with a 
finite spin polarization. In analogy to the edge states in 
the quantum spin-Hall effect, 15 we also consider the pos- 
sibility of edge states localized at the interface between 
regions with and without SOC. 

This paper is organized as follows. In Sec. II we intro- 
duce the model and the transfer matrix formalism used 
in the rest of the paper. In Sec. Ill we discuss the scat- 
tering problem at a single interface and the spectrum of 
edge states. In Sec. IV we address the case of a dou- 
ble interface — a SOC barrier — and the final Sec. V is 
devoted to the discussion of results and conclusions. 



II. MODEL AND FORMALISM 

We consider a clean graphene sheet in the xy- 
plane with SOCs 15,16 ' 21 ' 32,33 inhomogeneous along the 
x-direction. We shall restrict ourselves to a single- 
particle picture and neglect electron-electron interaction 
effects. The length scale over which the SOCs vary is as- 
sumed to be much larger than graphenc's lattice constant 
(a = 0.246 nm) but much smaller than the typical Fermi 



wavelength of quasiparticles Af. Since close to the Dirac 
points Ap ~ at low energy E this approximation 

is justified. This assumption ensures that we can use the 
continuum DW description, in which the two valleys are 
not coupled. Yet close to a Dirac point we can approxi- 
mate the variation of SOCs as a sharp change. Focusing 
on a single valley, the single-particle Hamiltonian reads 



H = v F a ■ p + Hsa, 
X(x) 

Hso = -rr 1 (<x X s), + A(x)<T z S 2 



(1) 
(2) 



where v-p w I0 6 m/s is the Fermi velocity in graphene. 
In the following we set h = v-p = 1 . The vector of Pauli 
matrices er = (a x ,a y ) [resp. s = (s x ,s y )] acts in sublat- 
tice space [resp. spin space]. The term Jiso contains the 
extrinsic or Rashba SOC of strength A and the intrinsic 
SOC of strength A. While experimentally the Rashba 
SOC can be enhanced by appropriate optimization of 
the substrate up to values of the order of 14 meV, 25 the 
intrinsic SOC seems at least two orders of magnitude 
smaller. Yet, the limit of large intrinsic SOC is of con- 
siderable interest since in this regime graphene becomes 
a topological insulator. 15 Thus in this paper we shall 
not restrict ourself to the experimentally relevant regime 
A ^> A but consider also the complementary regime. 
The wave function 'J is expressed as 

where the superscript T denotes transposition. Spectrum 
and eigenspinors of the Hamiltonian (I) with uniform 
SOCs are briefly reviewed in Appendix A. The spec- 
trum consists of four branches J5 aj£ (k) labelled by the 
two quantum numbers e = ±1 and a = ±1. Here, the 
first distinguishes particle and hole branches, the second 
gives the sign of the expectation value of the spin pro- 
jection along the in-plane direction perpendicular to the 
propagation direction k. The spectrum strongly depends 
on the ratio 



A 



(3) 



For rj > 1/2 a gap separates particle and hole branches. 
The gap closes at 77 = 1/2 and for r\ < 1/2 one particle 
branch and one hole branch are degenerate at k = (see 
Fig. 8 in App. A). 

We now briefly summarize the transfer matrix ap- 
proach employed in this paper to solve the DW scat- 
tering problem. 35-38 We assume translational invariancc 
in the y-direction, thus the scattering problem for the 
Hamiltonian (1) reduces to an effectively one-dimensional 
(ID) one. The wave function factorizes as ty(x,y) = 
e lkyV xi x )i where k y is the conserved y-component of the 
momentum, which parameterizes the eigenfunctions of 
the Hamiltonian of given energy E. 

For simplicity we consider pieccwise constant profiles 
of SOCs, and solve the DW equation in each region of 
constant couplings. Then we introduce the x-dependent 
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Fig. 1. Conservation of k y implies that 
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Figure 2: (Color online) Refraction angles as function of the 
incidence angle for fixed energy and fixed SOCs. Panel (a): 
E = 5, A = 0.5, A = 2; panel (b): E = 5, A = 2, A = 0.5. 



4x4 matrix f2(x), whose columns are given by the com- 
ponents of the independent, normalized cigenspinor of 
the ID DW Hamiltonian at fixed energy. 39 Due to the 
continuity of the wave function at each interface between 
regions of different SOC, the wave function on the left 
of the interface can be expressed in terms of the wave 
function on the right via the transfer matrix 



M = [n{x-)} 1 n( x +), 



(4) 



where xq is the position of the interface and x J = x<j ± 5 
with infinitesimal positive S. The condition det M. = 1 
guarantees conservation of the probability current across 
the interface. The generalization to the case of a sequence 
of N interfaces at positions Xi, i = 1, . . . , N, is straight- 
forward since the transfer matrices relative to individual 
interfaces combine via matrix multiplication: 



N 



m = H[n(x;)]- n(aH 



(5) 



From the transfer matrix it is straightforward to deter- 
mine transmission and reflection matrices, which encode 
all the relevant information on the scattering properties. 



III. THE N-SO INTERFACE 

First we concentrate on the elastic scattering problem 
at the interface separating a normal region N (x < 0), 
where SOCs vanish, and a SO region (x > 0), where 
SOCs are finite and uniform. 

We consider a quasiparticle of energy E, with E as- 
sumed positive for definiteness and outside the gap pos- 
sibly opened by SOCs. This quasiparticle incident from 
the normal region is characterized by the y-component 
of the momentum, or equivalently, the incidence angle 4> 
measured with respect to the normal at the interface, see 



„so 



E sin <j) = 
E coscj) 
E a cos £ Q 



E a sin £ Q = k. 



so 



(6a) 
(6b) 
(6c) 



where a = ±1 and E a = ^{E - A)(E + A - aXj. The 
refraction angles £ a are fixed by momentum conservation 
along the interface (6a) and read 



E 
E a 



■ sine 



(7) 



Figure 1 illustrates the refraction process at the N- 
SO interface. The incident wave function, assumed to 
have fixed spin projection in the z-direction, in the SO 
region splits in a superposition of eigenstates of the 
SOCs Hamiltonian corresponding to states in the differ- 
ent branches of the spectrum. These eigenstates prop- 
agate along two distinct directions characterized by the 
angles £ Q , whose difference depends on SOC and is an in- 
creasing function of the incidence angle, see Fig. 2. The 
angles £ a coincide only for normal incidence or for A = 0. 

Equation (7) implies that there exists a critical angle 
for each of the two modes given by 



~E 



(8) 



For <p larger than both critical angles <p a , the quasipar- 
ticle is fully reflected, since there are no available trans- 
mission channels in the SO region. For <f> in between the 
two critical angles the quasiparticle transmits only in one 
channel. 40 

After this qualitative discussion of the kinematics, we 
now present the exact solution of the scattering problem. 
In the N region x < a normalized scattering state of 
energy E > 0, incident from the left on the interface with 
incidence angle <j> and spin projection s is given by 

XN(aO=[*t,.|t>+<^.|4.>] ' 
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(9) 



where k x = k x (cf. Eq. 6b). Here the index s =t,l 
specifies the spin projection of the incoming quasipar- 
ticle with | f) and | J.) eigenstates of s z and dij is the 
Kronecker delta. The velocity Vp = cos0 is included to 
ensure proper normalization of the scattering state. The 
complex coefficients r s i s are reflection probability ampli- 
tudes for a quasiparticle with spin s to be reflected with 
spin s' . The associated matrix SIn^) reads 
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spin-down incident quasiparticle: 



Figure 3: (Color online) Angular dependence of the transmis- 
sion probabilities T+f (blue dashed line) and T-f (red solid 
line) at energy E — 2.5. The SOC are fixed as follows: (a) 
A = 0.1 and A = 0, (b) A = and A = 0.1, and (c) A = 0.5 
and A = 0.1. 



Similarly the wave function in the SO region (x > 0) can 
be expressed in general form as 



Xso(x) 



1 



\p+ip++(x) +r+ijj ++ (x)] 



+ -T=[M>-+(x)+r_$_ + (x)] (10) 

V v -+ 

where t± (resp. r±) are complex amplitudes for right- 
moving (resp. left-moving) states. The coefficient t a 
represents the transmission amplitude into mode a. The 
wave functions ip a + and the Fermi velocities v% + in the 
SO region are obtained from the expressions given in 
App. A with the replacement k x — > fcf^j where for no- 
tational simplicity the label SO will be understood. The 
wave functions ip a + are m turn obtained from ipa+ by 
replacing k xa — > —k xa . The matrix fisot^) then reads 



n so (a;) = 
e 2 



(11) 



1C 2 



\ie 



e 2 

. £± 
le 2 



e-'f— 
e"2" 
— \e~ 



9_ 

— ie~ 







o 



Af+e 













A/"_e ife 



N- 



\ 




,—ikx-x j 



where in the second matrix the normalization factors are 
defined as Af a = l/(2yjv a + sinh# Q ). 

According to Eq. (4) the transfer matrix for the sin- 
gle interface problem is given by the matrix product 
M = [f2 N (Cr)]~ 1 f2 so (0 + ). From At we obtain the trans- 
mission and the reflection probabilities for a spin-up or 
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(12) 
(13) 

(14) 
(15) 



where T a ((/>) = 9{4> a — <j))9(4> a + </>) with the Heav- 
iside step function. Here, T as is the probability for an 
incident quasiparticle with spin projection s to be trans- 
mitted in mode a in the SO region. Of course, probability 
current conservation enforces T +s + T_ s + R^ s + R^ s = 1. 

Figures 3 (a)-(c) show the angular dependence of the 
transmission probabilities for an incident spin-up quasi- 
particle into the (+) and (— ) modes of the SO region for 
different values of the SOCs. Panel (a) refers to the case 
of vanishing intrinsic SOC (A = 0). Here the (+) and 
the (— ) energy bands are separated by a SOC- induced 
splitting AE'cxt = A. Therefore at fixed energy the two 
propagating modes in the SO region have two different 
momenta, which gives rise to the two different critical 
angles (cf. Eq. (8) with A = 0). Panel (b) refers to the 
case A = 0, where the SOC opens a gap AE int = 2A 
between the particle- and the hole-branches, however the 
(+)/( — )-modes remain degenerate. Therefore at fixed 
energy these modes have the same momentum and, as a 
consequence, the same critical angles (cf. Eq. (8) for 
A = and A ^ 0). When both SOCs are finite - 
the situation illustrated in panel (c) — the transmission 
probabilities exhibit more structure. For incidence angles 
smaller than 4> + no particular differences with the cases 
of panels (a) and (b) are visible. When the (+) mode 
is closed, an increase (resp. decrease) of the (— ) mode 
transmission is observed for positive (resp. negative) an- 
gles, before the transmission drops to zero for incidence 
angles approaching 0_. The asymmetry between posi- 
tive and negative angles is reversed if the spin state of 
the incident quasiparticle is reversed. 

These symmetry properties can be rationalized by con- 
sidering the operator of mirror symmetry through the 
x-axes. 41 This consists of the transformation y —> —y 
and at the same time the inversion of the spin and the 
pseudo-spin states. It reads 



Sy = (a x ® Sy)Ry, 



(16) 



where R y transforms y — > —y. The operator S y com- 
mutes with the total Hamiltonian of the system [S y , Ho + 
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^so] = 0, therefore allows for a common basis of eigen- 
states. For the scattering states in the SO region (10) 
we have S y x+(Z+) = X+(M and S y X-(£-) = -%_(£_). 
Instead, it induces the following transformation on the 
scattering states (9) in the normal region: S y Xs{4>) = 
iX- s {~ 0)- By comparing the original scattering matrix 
with the ^-transformed one we find that 

T a , s (0)=T a ,_ s (-<£) (17) 

with a = ± and s =f, J,, which is indeed the symmetry 
observed in the plots. The asymmetry of the transmission 
coefficients occurs only when both SOCs are finite. 



A. Edge states at the interface 

In addition to scattering solutions of the DW equation, 
it is interesting to study the possibility that edge states 
exist at the N-SO interface, which propagate along the 
interface but decay exponentially away from it. The in- 
terest in these types of solutions is connected to the study 
of topological insulators. It has been shown — first by 
Kane and Mele 15 — that a zig-zag graphene nanoribbon 
with intrinsic SOC supports dissipationlcss edge states 
within the SOC gap. In fact, similar states are always 
expected to exist at the interface between a topologically 
trivial and a topologically non-trivial insulator. In our 
case, the latter is represented by graphene with intrinsic 
SOC. Of course SOC-frce graphene is not an insulator, 
however it is topologically trivial and edge state solutions 
do arise for \k y \ > \E\. When E is within the gap in the 
SO region the corresponding mode is evanescent along 
the x direction on both sides of the interface. Note that 
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Figure 4: (Color online). Energy dispersion of the edge state 
at the N-SO interface as a function of the momentum along 
the interface k y for different values of SOCs. Solution of the 
transcendental equation is allowed only for \k y \ > \E\ (white 
area). In all three cases shown n > 1/2: A = 1 and A = 0.4 
(lower-red line), A = 1.5 and A = 0.7 (middle-blue line), and 
A = 2 and A = 0.9 (upper-green line). 



the edge state we find is different from the one discussed 
in Refs. 15, 32 where zig-zag or hard-wall boundary con- 
ditions at the edge of the SOC region were imposed. 
The wave function on the N side then reads 

Xn(s)= (A\j)+B\t))e" x (18) 

with q = y/\k y \ 2 — E 2 . In the SO region the wave func- 
tion can be written as 
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with q a = ^k 2 - (E - A){E + A - aX). The continu- 
ity of the wave function at the N-SO interface leads to 
a linear system of equations for the amplitudes A to D. 
The matrix of coefficients must have a vanishing deter- 
minant for a non-trivial solution to exist. This condition 
provides a transcendental equation for the energy of pos- 
sible edge states, whose solutions are illustrated in Fig. 4 
for different values of the intrinsic and extrinsic SOCs. 
The condition \ k y \ > \E\ implies that solutions only exist 
outside the shadowed area. In addition, they are allowed 
only in the case SOCs open a gap in the energy spectrum, 
which occurs when 77 > 1/2 (see App. A and Eq. (3)). As 
can be seen in Fig. 4 the result is quite insensitive to the 
precise value of the extrinsic SOC. 

Edge states exist only for values of the momentum 
along the interface larger than the intrinsic SOC, i.e., 
k y > ky lln — A. The apparent breaking of time- reversal 
invariance (the dispersion is not even in k y ) is due to the 
fact that we are considering a single-valley theory. The 
full two- valley SOC Hamiltonian is invariant under time- 
reversal symmetry, that interchanges the valley quantum 
number. This invariance implies that solutions for neg- 
ative values of k y can be obtained by considering the 
Dirac-Weyl Hamiltonian relative to the other valley. The 
two counter-propagating edge states live then at opposite 
valleys and have opposite spin state and realize a peculiar 
ID electronic system. 

As mentioned in the Introduction, the intrinsic SOC 
is estimated to be much smaller than the extrinsic one, 
therefore in a realistic situation one would not expect the 
opening of a significant energy gap and the presence of 
edge states. It would be interesting to explore the pos- 
sibility to artificially enhance the intrinsic SOC, thereby 
realizing the condition for the occurrence of edge states. 



IV. THE N-SO-N INTERFACE 

The analysis of the scattering problem on a N-SO in- 
terface of the previous section can be straightforwardly 
generalized to the case of a double N-SO-N interface (SO 
barrier). Here the transmission matrix D is given by 
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Figure 5: (Color online). Panel (a): Angular plots for Tyt as 
a function of the injection angle for E — 2, A = 1 and A = 0. 
The three lines correspond to different distance between the 
interfaces: d — tt/2 (dashed black), d — n (dotted red), and 
d — 2n (solid blue). The spin-precession length is £so = f. 
When A = the transmission probability in the spin state 
opposed to the injected spin is always zero. Panel (b) and 
(c): angular plots of T-f-f (solid-blue) and (dashed red) as 
a function of the injection angle for E — 2, A = 1 and A = 0. 
The distance between the two interfaces is d — n in panel 
(a) and d — 2n in panel (b). The spin-precession length is 
£so = 2tt. 



Eq. (5) with N = 2. The transmission and the reflection 
probabilities in the case of a spin-up or -down incident 
quasiparticle read 
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(19) 
(20) 

(21) 

(22) 



In this case there is an additional parameter which con- 
trols the scattering properties of the structure, namely 
the width d of the SO region. In order to compare this 
length to a characteristic length scale of the system, we 
introduce the spin-precession length defined as 



= 2tt 



hv-p 

A + 2A 



(23) 



The intrinsic SOC alone cannot induce a spin preces- 
sion on the carriers injected into the SO barrier — an 
injected spin state, say up, is obviously never converted 
into a spin-down state. Figure 5(a) shows the angular 



Figure 6: (Color online). Angular plot of Ttt (solid-blue) 
and T±t (dashed-red) as a function of the injection angle for 
E = 2, A = 1 and (a) A = A/4, (b) A = A/2, and A = A. The 
distance between the two interfaces is kept fixed to d = £so- 



dependence of the transmission in the case of injection of 
spin-up. The behavior of the transmission as a function 
of the injection angle depends sensitively on the width d 
compared to the spin-precession length. For small width 
d < £so (dashed line) the transmission is a smooth de- 
creasing function of <p and stays finite also for <f> larger 
than the critical angle. In the case d > iso (dotted- and 
solid-lines) instead the transmission probability exhibits 
a resonant behavior and drops to zero as soon as the 
injection angle equals the critical angle. 

When only the extrinsic SOC is finite, the transmis- 
sion behavior changes drastically. Two different critical 
angles appear — the biggest coincides usually with 7r/2. 
The extrinsic SOC induces spin precession because of the 
coupling between the pseudo- and the real-spin. This is 
illustrated in Fig. 5(b)-(c). In Panel (b) we consider the 
case of spin-up injection with d = £so/2- At normal inci- 
dence the transmission is entirely in the spin-down chan- 
nel (dashed line). Moving away from normal incidence, 
the transmission in the spin-up channel (solid line) in- 
creases from zero and, after the first critical angle, the 
transmissions in spin-up and spin-down channels tend to 
coincide. In panel (c) the width of the barrier is set to 
d = £go- Here, the width of the SO region permits to an 
injected carrier at normal incidence to perform a com- 
plete precession of its spin state — the transmission is in 
the spin-up channel. For finite injection angles the spin- 
down transmission (dashed line) also becomes finite. For 
<t> < 0+ the transmission in the spin-up channel is almost 
fully suppressed while that in the spin-down channel is 
large. Finally, for <f> > (f> + the two transmission coeffi- 
cients do not show appreciable difference. 

In the case where both extrinsic and intrinsic SOC 
are finite, the transmission probability exhibits a richer 
structure. We focus again on the case of injection of 
spin-up quasiparticles. Moreover we fix the width of the 
SO region so that it is always equal to the spin-precession 
length d = £so ■ Fig- 6 illustrates the transmission proba- 
bilities T s t for three values of the ratio A/A = 1/4, 1/2, 1. 
Notice that from panel (a) to (c) the width of SO region 
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Figure 7: (Color online). Panel (a): total transmission T 
as a function of the injection angle for E = 2, d = 2-zr and 
several values of SOCs: A = 1 and A = (blue-solid line), 
A = and A = 0.5 (red-dotted line), A = 1 and A = A/4 
(yellow-dashed line), A = A/2 (orange-dashed-dotted line), 
and A = A (black-dotted-dotted-dashed line). Panel (b): z- 
component of the spin polarization Vz as a function on the 
injection angle for E = 2 and d = 2n and the following values 
of the SOCs: A = 1, A = and A = 0, A = 1 (same black- 
dashed line), A = 1 and A = A/4 (red-dotted), A = A/2 
(blue-dotted-dashed line), and A = A (green-solid line). 



decreases. 

The symmetry properties of the transmission func- 
tion can be rationalized by using the symmetry opera- 
tion (16). Proceeding in a similar manner as in the case 
of the single interface, for the SO barrier we find the 
following symmetry relations 

T s , s (ct>)=T S!S (-4>), (24a) 
T s _ a (<£) =T_ s , a (-<£), (24b) 

which are confirmed by the explicit calculations. 

So far we have considered the injection of a pure spin 
state — the injected carrier was either in a spin-up state 
or a spin-down state. Following Ref. 30 we now address 
the transmission of an unpolarized statistical mixture of 
spin-up and spin-down carriers. This will characterize 
the spin-filtering properties of the SO barrier. In the 
injection N region, an unpolarized statistical mixture of 
spins is defined by the density matrix 

pin = ^ IxtXxtl + glx+Xx-i-l. ( 25 ) 

where \ Xs ) = \s) ® \a) with \a) = (l/y/2)(l, e^) cor- 
responds to a pure spin state. When traveling through 



the SO region, the injected spin-unpolarizcd state is sub- 
jected to spin-precession. The density matrix in the out- 
put N region can be expressed in terms of the transmis- 
sion functions (19) as 

Pout = ir T |Ct)(Ctl + iniCi)(Cil, (26) 

where the coefficients T s = Tf s + T± a are the total trans- 
missions for fixed injection state. The spinor part is de- 
fined as 

•<■>-;&: (£)•">■ <27) 

where t s : s arc the transmission amplitudes for incoming 
(resp. outgoing) spin s (resp. s'). The output density 
matrix is used to define the total transmission 



and the expectation value of the z component of the spin- 
polarization 

v * = \ ( T n + T n - Tit - Tu) ■ (29) 

In Fig. 7 we report the total transmission (panel (a)) 
and the z-component of the spin-polarization (panel (b)) 
as a function of the injection angle for fixed energy and 
width of the SO region. We observe that for an un- 
polarized injected state the transmission probability is 
an even function of the injection angle T(4>) = T{—<f), 
Moreover, for injection angles larger than the first criti- 
cal angle <f> > <f> + , the transmission has an upper bound 
at T = 1/2. On the contrary V z is an odd function of 
the injection angle V z ((f>) = —V z {—(f>). It is zero when 
at least one SOC is zero. When both SOC parameters 
are finite V z is finite and reaches the largest values for 
<f> > <j>+- The maxima in this case increase as a function 
of the intrinsic SOC. 

To experimentally observe this polarization effect the 
measurement should not involve an average over the an- 
gle 4>, which, otherwise — due to the antisymmetry of V z 
— would wash out the effect. To achieve this, one could 
use, e.g., magnetic barriers, 37,42 which are known to act 
as wave vector filters. 



V. CONCLUSIONS 

In this paper we have studied the spin-resolved trans- 
mission through SO nanostructures in graphene, i.e., sys- 
tems where the strength of SOCs — both intrinsic and 
extrinsic — is spatially modulated. We have considered 
the case of an interface separating a normal region from a 
SO region, and a barrier geometry with a region of finite 
SOC sandwiched between two normal regions. We have 
shown that — because of the lift of spin degeneracy due 
to the SOCs — the scattering at the single interface gives 
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rise to spin-double refraction: a carrier injected from the 
normal region propagates into the SO region along two 
different directions as a superposition of the two avail- 
able channels. The transmission into each of the two 
channels depends sensitively on the injection angle and 
on the values of SOC parameters. In the case of a SO 
barrier, this result can be used to select preferential di- 
rections along which the spin polarization of an initially 
unpolarizcd carrier is strongly enhanced. 

We have also analyzed the edge states occurring in the 
single interface problem in an appropriate range of pa- 
rameters. These states exist when the SOCs open a gap 
in the energy spectrum and correspond to the gapless 
edge states supported by the boundary of topological in- 
sulators. 

A natural follow-up to this work would be the detailed 
analysis of transport properties of such SO nanostruc- 
tures. From our results for the transmission probabil- 
ities, spin-resolved conductance and noise could easily 
be calculated by means of the Landaucr-Biittiker formal- 
ism. Moreover we plan to study other geometries, as, e.g., 
nanostructures with a periodic modulation of SOCs. The 
effects of various types of impurities on the properties 
discussed here is yet another interesting issue to address. 

We hope that our work will stimulate further theo- 
retical and experimental investigations on spin transport 
properties in graphene nanostructures. 
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C e (k) = 



1 



2\/cosh 9 a 



(A2) 



(e -i*- £ Q / 2; ee ee a /2^ iaee ^/ 2 , iae^"^ 2 ), 



where T denotes transposition and 

sinh# Q = — — , 

k 

as k x + ik y 



(A3) 
(A4) 



with k = ^Jk% + k^. The spin operator components are 

expressed as Sj = \sj <g> cq. Their expectation values in 
the eigenstate <E> Qie read 

<Sx> 
(Sy) 



—ea sin <p 
2 cosh 6 a ' 
ea cos <b 



2 cosh t 
= 0, 



(A5a) 

(A5b) 
(A5c) 



which shows that the product ea coincides with the sign 
of the expectation value of the spin projection along the 
in plane direction perpendicular to the direction of propa- 
gation. For vanishing extrinsic SOC, the eigenstates $ a ,e 
reduce to linear combinations of eigenstates of S z . 

Similarly, the expectation value of the pseudo-spin op- 
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Appendix A: Graphene with uniform spin-orbit 
interactions. 



In this appendix we briefly review the basic proper- 
ties of DW fermions in graphene with homogeneous SO 
interactions. 21 The energy eigenstates are plane waves 
ip ~ $(k)c lk ' r with $ a four-component spinor and eigen- 
values given by («p = ft = 1) 




(Al) 



where a = ± and e = ±. The energy dispersion as a 
function of k x at fixed k y = is illustrated in Fig. 8 for 
several values of A and A. The index e = ± specifics the 
particle/hole branches of the spectrum. The eigenspinors 



(a) 




(c) 





-(b) 




-2 2 

Momentum 



-2 2 

Momentum 



Figure 8: Spectrum of the DW Hamiltonian with intrinsic 
and Rashba SOC as a function of k x for k y — . Panel (a): 
dashed lines refer to A = 0.5 and A = 0; solid and dotted lines 
refer to A = and A = 1. Panel (b): A = 0.4 and A = 1. 
Panel (c): A = 0.5 and A = 1. Panel (d): A = 0.8 and A = 1. 
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crator cr is given by 



ecoscp 
cosh 9 a 

esincf) 
cosh 6 n 



(A6a) 
(A6b) 



Since the SOCs in graphene do not depend on momen- 
tum, the velocity operator still coincides with the pseudo- 
spin operator: v = r = i[H, r] = cr. Thus the velocity 



expectation value in the state $ Qj£ is given by Eqs. (A6a- 
A6b). Alternatively, it can be obtained from the energy 
dispersion as 



V k £ Q 



fc 2 + (A - ^ 



(A7) 



The group velocity is then independent of the modulus of 
the wave vector if either the SOCs vanish or A = aX/2. 
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